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The multi-band Gutzwiller method, combined with calculations based on den- 
sity functional theory, is employed to study total energy curves of the ferro- 
magnetic ground state of Ni. A new method is presented which allows flow 
of charge between d and s, p type orbitals in an approximate way. Further it 
is emphasized that the missing repulsive contribution to the total energy at 
large magnetic moments can be estimated from an analysis of specific DFT 
calculations. 



1. INTRODUCTION 

Combining the multi-band Gutzwiller methodS and density functional 
theory (DFT) allows to study magnetic states of crystalline solids!'!. The 
presently used Gutzwiller DFT scheme employs a tight-binding model for 
the single particle energies derived from DFT calculations for non-magnetic 
crystal and uses certain chemical potentials which keep the partial charges 
n s , n p , and of s, p, and d orbitals constant for all values of the magnetic 
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moment /i. Calculations for ferromagnetic iron group metals, in particular 
for Ni, yield total energy versus magnetic moment curves Etot(fJ'), which are 
in good agreement with experiment, yet the values \x m in at minimum total 
energy are somewhat too large. 

In this paper we first present an alternate method which, in magnetic 
cases, allows charge flow between the various orbitals. We also argue that 
additional repulsive contributions to the total energy exist near and beyond 
the limit of strong ferromagnetism, which are not included in our scheme. 
These contributions can be estimated from an analysis of specific DFT cal- 
culations. When incorporated in total energy curves, we obtain [i m i n values 
very close to experiment. 



2. GROUND STATE ENERGY FOR NICKEL FROM THE 
MULTI-BAND GUTZWILLER METHOD 

2.1. General treatment 

We start from the Hamiltonian 

cr,cr f 

= H±+H at (1) 

Here, the quantities tf? are the matrix elements of the tight-binding Hamil- 
tonian Hi. For Ni we use a basis of 4s, 4p, and 3d orbitals, in total 2x9 = 18 
spin-orbitals a. The tf? values have been obtained from a least squares fit 
to the energy bands of a DFT calculation with the LAPW-WIEN code for 
non-magnetic Ni, using the local density approximation (LDA)i. In the 
fits we have preferably used energies of k points at high symmetry points 
and along high symmetry lines. This was done in such a way that the full 
information on the symmetry of all states involved in the fit was taken into 
account!!. The procedure allows to incorporate information on the orbital 
character of the states. 

The atomic Hamiltonians 

tt( 1 ) _ rr<xi...o-4 *f p.] p p 

<7l...(74 

= E E Hr)<r| (2) 

r 

are assumed to be the same for each lattice site i and include all on-site 
Coulomb interactions within the 3d shell. We use the spherical atom approx- 
imation so that all interactions are determined by the three Slater-Condon 
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integrals F°, F , F 2 or, equivalently, by the three Racah parameters A, B, 
C as outlined, e.g. in the textbook of Sugano et al i. 

In the subspace of the 3d shell, the states |T) represent all 2 2 ' 5 — (2-5+1) 
multi-electron eigenstates with eigenenergies Ey- 
The expectation value of H is given byEI 



**7 

cf.cf' 



icr iT 



$o) • (3) 



®0 \H e ff 

Here, 

= n ia = ($ |n«r| $o) (4) 

is the density in the spin-orbital a (we assume a monoatomic cubic lattice 
and skip all labels i, j, if possible) and e CT = = tff. The expectation value 
(jH) has to be minimized with respect to the single particle product wave 
function |$q) an d t ne occupancies mr of the atomic multi-electron states 
|r). Both determine the 'hopping reduction factors' q a via the relations 

(5) 

where the coefficients Sf r/ depend on the densities n a > and the eigenstates 
|T) of H^t ■ l^o) is the ground state of H e ff, but it has to be determined in 
a self-consistent way, since £f e // depends on |<3?o) via the densities n a . 

For the optimum set of variational parameters, the effective single parti- 
cle Hamiltonian leads to quasi-particle energy bands of a Fermi liquid which 
determine, e.g., the shape of the Fermi surface or may be compared to ex- 
perimental energy bands, obtained from angular resolved photoemission. 

The wave function |$o) can be chosen to yield magnetic ground states by 
incorporating orbital exchange splittings Ai2# and Ae g to produce majority 
and minority bands. Then ^ff'*) is a function of the magnetic moment \x. 
In fact, for a given value of the magnetic moment Et t(At2 g , Ae g ) shows 
a minimum, when At2 g /Ae g ~ 2-3; i.e. the exchange splitting is strongly 
anisotropic. 



2.2. The method of chemical potentials 



The unlimited variation of the spin-orbital charge densities n a leads 
to a charge flow from 3d to 4s and 4p states, since the Hamiltonian of eq. 
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(jTj) includes only atomic 3d electron-electron interactions. The charge flow 
can be suppressed by appropriate chemical potentials which keep the values 
of the charge densities fixed to specific values n^ FT , n p FT , and thus also 
n^ FT . These values are obtained from the ground state of Hi, the tight- 
binding model, which is derived from the DFT results for non-magnetic 
Ni. The choice is justified by the observation that, in general, charge den- 
sities obtained from DFT calculations show very good, albeit not perfect, 
agreement with experimental data. Typical total energy curves Et t([J>) for 
Ni are shown in Fig. [l]a. The first calculation has been carried out us- 
ing A = 10 eV, C = 0.4 eV and C/B = 4.5. /j, min « 0.68 is found, and 
Econd = E tot (fi=0) - Etot((irnin) ~ 40meV. Note that the ratio C/B is typi- 
cal for 3d ions and that the exchange interaction parameter C = 0.4 eV is an 
estimate derived from free Ni ionstl. The Racah parameter A corresponds 
to the Hubbard model parameter • U ; its value is adjusted to produce the 
experimental <i-band width of Ni Bu. 

The second calculation employs a significantly larger C = 1.0 eV, again 
with C/B = 4.5 and A = 10 eV. Here we obtain \i m in = 0.77 and E con d 
84meV. 

Our value for the Racah parameter A = 10 eV is much larger than 
suggested in other works (see, e.g., Refs. |7]]8|). This discrepancy may have 
to do with the neglect of hybridization terms between s, p and d orbitalsQ. If 
we switch off these terms in our theory we find A ~ 3 eV which is comparable 
to the results in Ref. pL 



2.3. Approximate incorporation of other contributions to the 
electron-electron interaction: the charge flow method 

The charge flow can be allowed if additional interaction terms involving 
the s and p type charge densities are included in the total energy. A simple 
term is of the form 

E sp = A sp<d (n s + n p ) n d + 7^ A s P (n s + n p ) 2 (6) 

If we add E sp to the Gutzwiller ground state energy for paramagnetic 
nickel, now unrestricted in the charge densities, we can choose A sp ^ and 
A sp in such a way as to yield the minimum total energy for the densities 
n^ FT and n p FT . Typical values of A sp ^ and A sp are of order of the Racah 
parameter A. Introducing magnetic states \&o(At2 g , Ae 9 )}, the changes in s- 
d and p-d hybridization due to the exchange splittings At2 g , Ae 9 of majority 
and minority bands leads to increasing n s , n p values (and decreasing n^). 
These changes cause E sp to increase and compensate for the reduction of the 
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Fig. 1. Gutzwiller DFT total energy curves versus magnetic moment [i 
for ferromagnetic Ni. a) shows calculations using the chemical potential 
scheme (exchange interactions C = 0.4 eV and C = 1.0 eV, respectively) 
and using the new charge flow method (C = 0.4 eV). Also shown is the 
LDA total energy curve, b) exhibits the repulsive part of the LDA total 
energy curve E^q A , obtained by keeping only the charge-dependent part of 
the exchange correlation potential. Further shown is the one-particle band- 
splitting contribution at small moments fi and the residual repulsive 
energy curve E^ A . c) shows the three Gutzwiller DFT total energy curves 
of a), corrected by either E^ A (for the C = 0.4 eV calculations) or by Eqq A 
(for the C = 1.0 eV calculation). 



T. Ohm, S. Weiser, R. Umstatter, W. Weber, and J. Biinemann 



d-d interaction energy (which in first approximation is \ A (n^ 2 ). A typical 
total energy curve (using A = 10 eV, C = 0.4 eV and C/B = 4.5) is also 
shown in fig. ||a. The curve is quite similar to the corresponding Etotin) 
curve using the chemical potential method, yet the value Umin f« 0.625 \ib is 
somewhat smaller than before. 



3. ANALYSIS OF DFT TOTAL ENERGY CALCULATIONS 

Using spin-dependent DFT, again within the LDA, we reproduce the 
well known E^ A {p) curve with its minimum at /W^iri = 0.62 /ig (also see 
fig. |l|a) . If we put equal to zero the spin-dependent part of the LDA exchange 
correlation functional, we obtain a repulsive curve Ej^, A Charge0nly (fi) (fig. 
|l]b). This curve E^ A (v) has its minimum at /i = and coincides there 
with E^ f A (fi = 0), as expected. The function E^Q A (fi) can be interpolated 
very well by the polynomial 

E L C D A ^) - E t L o ? A (0) = a 2 f + + <V, (7) 

with a 2 = 0.090 eV//i%, a 4 = 0.054 eV/n%, a 6 = 0.244 eV/(i%. 

The term a 2 y? dominates Eqq in the range values < \x < 0.4; i.e. 
close to the limit of strong ferromagnetism, where the majority d bands have 
been filled completely. The term a 2 (j? appears to arise completely from the 
splitting of the majority and minority bands. This can be shown by evalu- 
ating the sum over the LDA valence band energies. For the non-magnetic 
case, the tight-binding model of H\ agrees very well with the LDA energy 



bands i.e., E bs = H x <f> K x >), with 

being the ground state wave 

function of H\, gives the sum of all occupied LDA valence band energies. 
For small values of the exchange splittings At 2g , Ae g (and thus small values 
of the magnetic moment fi) the majority and minority LDA ci-bands shift 
rigidly, which is very well reproduced by corresponding ^^'(At 2 g, Ae g 
states. Thus, for small [a, the expression 



EM = {*a\*hg, Ae 9 ) H x &o>(At 2g , Ae g ) 

= 4 1) +a£V + ~ (8) 

should be compared to Ef J a f A ' CO (fi). Indeed, we find a^p ~ a 2 . However, 
the term Eb s (fi) i.e. the increase in total energy due to the splitting of 
minority and majority bands is included in the Gutzwiller treatment eq.(|3|). 
Therefore, the repulsive term 

E rep A (v) = + a 6 /i 6 (9) 
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Fig. 2. Minority quasi particle bands for ferromagnetic Ni (/x = 0.55 hb), 
obtained from the Gutzwiller DFT calculations using either C = 0.4 eV (full 
lines) or C = 1.0 eV (dashed lines). The two calculations differ only on a 
scale of a few meV. 
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Fig. 3. Majority quasi particle bands for ferromagnetic Ni (p, = 0.55 /j,b) 
from the C = 0.4 eV (full lines) and the C = 1.0 eV calculations (dashed 
lines). The two calculations differ on a scale of 10-100 meV. 
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is interpreted to represent the missing repulsive part of our Gutzwiller treat- 
ment. We remind our readers that the Hamiltonian of eq. ([l]) is still a very 
simplified model. It uses a very limited local basis, assumes the orbitals of 
the neighboring atoms to be orthogonal, only valence electrons are consid- 
ered and only part of the electron-electron interaction is included. On the 
other hand, density functional theory treatments have a rather complete ba- 
sis, incorporate all core states and all electron-electron interactions, albeit 
on the basis of an effective single particle theory. 

If we add E^® A to the total energy curves, obtained either using the 
chemical potential method or the charge flow method, we find very good 
agreement of fJ, m in with the experimental spin-only value of 0.55 \xb (see fig. 

Be). 

Alternatively, if we interprete the full curve Eqq A (/j,) as the repulsive 
part missing in the Gutzwiller treatment, the curves using C = 0.4 eV do 
not give any minimum for magnetic states. Only when we enhance C to 
C = 1 eV, a reasonable value of E con d is found, again for /x m m 

close to the 

experimental value. 

4. QUASI-PARTICLE ENERGY BANDS 

The change of C from 0.4 eV to 1.0 eV may not only produce quite 
different total energy curves, but might also result in drastic changes of the 
quasi particle energy bands. Therefore in fig. || and fig. |3| the corresponding 
bands for C = 0.4 eV, and C = 1.0 eV are compared, both calculated at 
H = 0.55 fiB and after minimizing the total energy with respect to the ratio 
At2 g /Ae g (which is rather similar for the two cases). We find few differences 
of the quasi particle bands, very little differences concerning the minority 
bands (see fig. |2|) and somewhat bigger ones for the majority bands (see fig. 
||). This result is not surprising, since even for C = 1.0 eV, we are still in the 
limit C/A <C 1, where the q a - values which enter the effective single particle 
Hamiltonian are mainly determined by A. The quasi-particle bands are also 
influenced by the exchange splittings Ai2 9 and Ae 9 . Again, since we find 
similar ratios Ai2 9 /Ae g for the two cases, the anisotropy of the exchange 
splitting is similar. 

5. CONCLUDING REMARKS 

We have presented an alternate scheme for carrying out Gutzwiller cal- 
culations which allows charge flow. This scheme uses simple approximations 
to treat in the total energy expression those electron-electron interaction 
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terms, which are not included in the Gutzwiller Hamiltonian. The interac- 
tion parameters are determined from paramagnetic calculations. This new 
scheme appears to be somewhat superior to the previously introduced chem- 
ical potential method, where all charges are kept constant. In both methods, 
the magnetic moment values at minimum total energy are somewhat larger 
than the corresponding experimental spin-only moment. 

We have tried to estimate the missing repulsive contribution to the 
total energy from an analysis of specific DFT calculations, where the spin- 
dependent part of the exchange-correlation potential has been switched off. 
When this repulsive term is corrected for the band splitting contributions 
already included in the Gutzwiller method, the resulting total energy curves 
yield values of the magnetic moments which agree very well with experiment. 

If we assume that the repulsive part of the LDA total energy curves 
remains uncorrected, i.e. if the band-splitting contribution is kept in the 
repulsive term, calculations using a bigger interaction parameter (7=1 eV 
lead to very satisfactory results concerning the magnetic moment. We then 
compare the quasi particle energy bands based on the two calculations with 
C values of 0.4 eV and 1.0 eV, respectively. Both calculations yield the same 
(experimental) value of the magnetic moment. We find very little differ- 
ences in the energies of minority bands and only moderate differences for 
the majority bands. 

These results justify a simplified Gutzwiller DFT treatment, where the 
magnetic moment is fixed to the experimental spin-only value and Etot is 
minimized with respect to the ratio At2 g /Ae g of the orbital exchange fields. 
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